Fidelity and dynamics near Dirac points in two-dimensional systems 
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We study quantum fidelity and dynamics near quantum critical points (QCPs) of the two- 
dimensional (2-D) Dirac Hamiltonian of graphene (the gapped to gapless transition induced by 
a mass term) and the 2-D BHZ Hamiltonian of HgTe/CdTe quantum wells (the topological to 
trivial insulator transition). For the two-dimensional Dirac Hamiltonian, we encounter marginal 
behaviour of the ground state fidelity near the Dirac point, which is displayed in the absence of a 
sharp dip in the ground state fidelity (or equivalently the weak logarithmic divergence of the fidelity 
susceptibility) . There is also a logarithmic correction to the proposed scaling of fidelity in the ther- 
modynamic limit. We then study the dynamics of the edge states of the 2-D BHZ Hamiltonian in 
a ribbon geometry following a sudden quench to the QCP. The effective edge state Hamiltonian is 
a collection of decoupled two-level systems which get coupled to bulk states following the quench. 
We notice a pronounced collapse and revival of the Lochschmidt echo for low-energy edge states 
illustrating the oscillation of the state between the two edges. We also observe a similar collapse 
and revival in the spin Hall current carried by these edge states, leading to a persistence of its 
time-averaged value. 
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The one-dimensional (1-D) and two-dimensional (2-D) 
Dirac Hamiltonians (DHs) have found a wide range of 
applications in quantum condensed matter systems in re- 
cent years PUS] . The low-energy physics of graphene [4] , 
and of the bulk states in 2-D topological insulators is 
described by a 2-D DH [3J |B] , and the 1-D edge states 
existing in a 2-D topological insulator (TI) system are 
described by an effective 1-D DH [3 . Quantum critical 
points (QCPs) [THS] of graphene (the gapless to gapped 
transition induced by a mass term) as well as the 2-D 
BHZ Hamiltonian (the topological insulator to trivial in- 
sulator transition) are 2-D Dirac points (DPs) with lin- 
ear dispersions. Experimentally, a gap can be opened in 
the otherwise gapless linear band structure of graphene 
through several methods e.g., by the application of an ex- 
ternal transverse electric field to the graphene sheet [TU] ; 
this is mimicked by adding a mass term to the DH which 
vanishes at the gapless QCP. The experimental prospect 
of real-time tuning of parameters controlling these tran- 
sitions in optical [H] [12] and photonic lattices [HI [14] as 
well as by exploiting Floquet dynamics [PHUT] , has made 
the study of these Hamiltonains from the viewpoint of 
quantum dynamics timely and necessary. 

In this paper, we study the ground state quantum 
fidelity [15] [15] of a 2-D DH and show the difference 
in results as compared to the 1-D case. We encounter 
marginality in the behavior of the ground state fidelity 
near the DP which shows up in a prominent way both in 
the scaling of the fidelity susceptibility and the scaling of 
the fidelity in the thermodynamic limit. The fidelity sus- 
ceptibility is proportional to the density of defects gener- 
ated by a sudden quench across a QCP [501 125 • We then 
focus on the dynamics of edge states of the BHZ Hamil- 
tonian when the system is suddenly quenched from the 
TI phase to either the QCP or the trivial insulator (TrI) 



phase. The quench couples the two-level subspace of the 
edge states to a multi-level environment of bulk states. 
We observe a strong oscillation of the low-energy edge 
states between the two edges of the system and the per- 
sistence of the spin Hall current (SHC) when the system 
is quenched to the QCP, which we attribute directly to 
the linear low-energy dispersion at this point (also found 
in other models describing 2-D TIs). 

The ground state quantum fidelity [T5] (F), which 
measures the overlap between many-body ground states 
at slightly different values of a parameter to of the 
Hamiltonian, is an important tool for detecting quan- 
tum phase transitions [22 - 25]: one can use the expansion 
F = \(ip(m + 5)\tl){m))\ 2 = 1 - 6 2 L d XF(m) + where 
5(— > 0) is a small change in the parameter to and L is the 
linear dimension of the c£-dimensional system, and xf(w) 
is defined to be the fidelity susceptibility density at to. 
The quantum fidelity generally displays a marked drop 
as to approaches the value m c at which the Hamiltonian 
has a QCP. Likewise, xf usually diverges in a universal 
power-law fashion as \m — m c \ 2 ~ vd where v is correlation 
length exponent associated with the QCP 22- 24 , I26]I27|. 
We explore the marginal case with vd = 2, using the ex- 
ample of the 2-D DH with a mass term, while the 1-D 
case has been found to satisfy the proposed scaling rela- 
tion with v = d — 1 28 . 

This Hamiltonian can be written in a two-state basis 
as follows: 

where the parameter m is the Dirac mass. This produces 
two bands with dispersions E± = ±^Jk 2 + k 2 + m 2 , and 

displays a QCP (DP) at m — (with exponent v = 1) 
where the gap between the two bands vanishes. 
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We obtain the fidelity F using the two-level nature of 
the Hamiltonian ([I]); for small values of S, we find 



= CXph T67 ( 
L 2 m 2 
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+ ln 
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L 2 7 



4tt 2 
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where k max is an upper cutoff on the k modes. Fidelity 
drops to zero when k max — > oo which shows that the high 
energy modes contribute significantly, which is charac- 
teristic of a marginal case (see Fig. [T]) . The logarithmic 
scaling of the fidelity susceptibility with m and L also 
points to the marginality of the situation. 
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FIG. 1: (a) Fidelity as given by (2}, calculated for 5 = 0.001 
and L X:V — 10,000. The fidelity drops to zero at the QCP 
(m = 0) but a sharp dip is absent, (b) The fidelity suscepti- 
bility density xf has a peak at the QCP. 

The most prominent role of marginality manifests in 
the fidelity studied in the thermodynamic limit (L — > oo 
and small but finite S), where one cannot use the fidelity 
susceptibility approach [29] . It has been conjectured that 
in this limit, lnF ~ -5 vd L d at the QCP; this scaling 
has been verified for the 1-D DH [2S]. In the 2D case, 
we obtain \nF = — L 2 6 2 In §/(8tt). We therefore have a 
marginal logarithmic correction to the expected 5 2 scal- 
ing of the fidelity even in the thermodynamic limit. 

We now turn to the 2-D BHZ Hamiltonian describing 
the low-energy electrons in Hg-Te/Cd-Te quantum well 
TIs, which displays helical edge state solutions that exist 
within the bulk bandgap in the TI phase [3] |6]. The 
bulk states undergo a QPT with the low-energy modes 
satisfying a 2-D DH with a linear dispersion at the QCP 
(2-D DP). Additionally, there are chiral edge states with 
linear dispersion in the TI phase which are described by 
effective 1-D DH [3J. We study the decoherence of these 
edge states using the Lochschmidt echo (LE) [30J when 
the 2-D Hamiltonian is quenched from the TI phase to 
the QCP (or to the TrI phase). 

The 4x4 Hamiltonian comprising of two 2x2 blocks 
(for opposite electron spins) is given by 
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FIG. 2: Spectrum of the BHZ hamiltonian in a ribbon geome- 
try with ribbon width L — 200 nm. The parameters used are 
A = 364.5 meV/nm, B = -686 meV/nm 2 and C = D = 0. 
The values of m used are —10 meV (a), (b) and +10 meV 
(c). There is a small gap of 0(A/L) in the spectrum near 
E = at m = because of the finite width of the ribbon. 
There also exists an exponentially small gap between the edge 
state bands in the TI phase. Note the (almost) equal spacing 
of energy levels for k x = at m = 0, characterized by the 
solutions of a Dirac particle in a 1-D box. 



[m — B{k 2 + k 2 )]a z . Here, A,B,C,D and m are de- 
termined by the thickness of the quantum well and the 
material parameters; the parameter m controls the phase 
of the system and changes sign relative to B when the 
system crosses from the TI phase (where edge states 
are present) to the TrI phase (with no edge states) via 
a DP at m I B = 0. This Hamiltonian differs from 
the previously discussed Hamiltonian by the term 
—B(k 2 + k 2 )a z , which allows for the existence of edge 
states which cross the bulk bandgap. Although the re- 
sults presented here are valid in generic situations, we 
shall set D = for simplicity, which also ensures an 
electron-hole symmetric spectrum. We consider a ribbon 
geometry extending from —L/2 to L/2 in the y direction 
(with the wavefunction vanishing at the edges) and apply 
periodic boundary conditions in the a;-direction |31H33j . 

Analyzing Hamiltonian ([3]), with the given boundary 
condition, one can show that in the TI phase (m/B > 0), 
there are two types of eigenstates of the Hamiltonian, 
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i.e., edge states (localized towards the edges and decay- 
ing over a length 1/A(fc x ,m)) and bulk states (spreading 
across the whole ribbon). The spectrum, which is sym- 
metric in ±k x and ±E, is displayed in Fig. [2j The edge 
states in the TI phase exist for (fc^l < ko, fco depends on 
m and L [52] . 

We perform a sudden quench of the parameter m going 
from m/B>0tom/B<0 and look at the subsequent 
evolution of an edge state and its spin current. The edge 
states, which originally formed a qubit-like two level sys- 
tem with an effective Hamiltonian H ec i ge rs Ak x a z at 
each k x [3] now get coupled to several bulk modes and 
subsequently decohere. Following a sudden quench, the 
evolution of an edge state is given by 



\i>edge{k x ,t)) 



E 

n=—oc 



(lpn{k x )\lpedge(k x ))e | lj) n (k x )) , 

(4) 

where \ip e dge(k x )) is an edge eigenstate of the Hamil- 
tonian at the initial value of m (= mi) and \tp n (k x )) 
are the eigenstates of the Hamiltonian at the final 
value m%. The index n runs from — oo to oo exclud- 
ing n = and denotes the -ve and +ve energy bulk 
modes, respectively. Since all the modes are plane waves 
along the x direction, different k x modes do not cou- 
ple to each other. To study the dynamics of a single 
edge state and quantify its decay, we calculate the LE 
C{t) = KVWgele^" 11 ^-^™ 2 ) 4 !^)! 2 , which using 
Eq.Q can be put in the form 



oo 

E 

n=— cx 



\(ipn(k x )\ip ed ge(k x ))\ 2 e~ 



iE n t 



(5) 

In general, the LE defined above initially drops rapidly 
with time and turns into a rapidly oscillating noisy func- 
tion of small amplitude (Fig. [3]ja)), indicating that the 
edge state decoheres significantly. However, there is a 
striking difference when one looks at the evolution of a 
low-energy (k x « k ) edge state following a quench to 
the QCP at m = 0; the LE of the edge state shows a 
pronounced collapse and nearly complete revival for sev- 
eral oscillation cycles (Fig. [3](b)) . This is a consequence 
of the nearly equal spacing of the first few energy levels 
at low k x near m — (arising due to confinement of the 
linearly dispersing particles in a ribbon geometry) where 
the overlap with the edge state is the most significant (see 
Fig. [2). We then have E n « sign(n) [E g /2 + (|n| - l)AE] 
for all significant terms in Eq.Q, where E g ~ 1/L is the 
bulk bandgap. The summations over n > and n < 0, 
then represent Fourier series of a peroidic function with 
period r = 2tt/AE making the LE a periodic function 
with this period. Since C(t = 0) = 1, the LE shows 
a near-complete revival at t = nr. Eventually, after 
several oscillations, the slight non- uniformity in spacing 
becomes significant and the revival of the LE weakens. 
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FIG. 3: Lochschmidt Echo for various edge states and 
quenches. The system parameters are A = 364.5 meV/nm, 
B = -686 meV/nm 2 , C = D = and L = 400 nm. LE for an 
edge state with: (a) k x = 0.01 meV/nm and m = — 10 meV 
after quenching to m — +10 meV. There is no significant 
revival of the edge state, (b) k x — 0.001 meV/nm and 
m = —10 meV after quenching to the QCP at m = 0. There 
is a pronounced collapse to and a nearly full recovery of the 
LE for several cycles, (c) k x = and m = — 10 meV after 
quenching to the QCP. There is a doubling of the frequency 
of oscillation as compared to the previous case, as this edge 
state exists on both edges. 



Since AE ~ the period of this revival scales as L. 
For small quench amplitudes (\m\ << A/L), the edge 
state does not decay significantly. 

Interestingly, we find that the edge state travels from 
one edge to the other and back, existing on opposite 
edges at the points of maxima and minima of the LE 
(Fig. Ha)). This effect is due to the finite width of the 
ribbon and will not be seen in an infinite system. Since, 
for a significantly large system, the edge states at oppo- 
site edges do not overlap, the LE drops to zero when the 
edge state reaches the opposite edge, and revives again 
when it comes back. For very low values of momentum 
k x — > 0, the edge state exists with peaks on both edges 
of a finite ribbon |31| . Hence, when a peak on a given 
edge travels to the opposite edge, the peak on the oppo- 
site edge also travels simultaneously to the given edge, 
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FIG. 4: (a) The probability density p — ip^ip of an edge state 
(k x = 0.001 meV/nm) following a sudden quench from m = 
— 10 meV to m = shown at t = (solid) and t = r/2 
(dashed). The edge state travels between the two edges, (b) 
The probabilty current density of the state in the x direction 
(J x ) at t = (solid) and t = r/2 (dashed). The same state 
carries currents of opposite direction on opposite edges. The 
system parameters are the same as Fig. j3b. 



resulting in a maximum of the LE at (2n + l)r/2 instead 
of a minimum, and hence a doubling of the frequency of 
oscillation of the LE. The LE is now minimum at times 
when the state is concenterated near the middle of the 
ribbon (Fig. go)). 
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FIG. 5: The net probabilty current carried by the pair of edge 
states with E = 3.65 meV (near the —L/2 edge over a length 
1/A ~ 37 m) following a quench from m = — 10 me7 to m = 
0. The current collapses and revives like the LE, and is always 
greater than zero, leading to a persistence of its time-averaged 
value. The system parameters are A — 364.5 meV/nm, B = 
-686 meV/nm 2 , C = D = and L = 200 nm. 



The probability current carried near the edge in the x 
direction by the edge states over their decay length 1/A 
is proportional to the net SHC carried by the state and 
its time-reversed conjugate in the opposite spin sector. 



It can be calculated using the continuity equation for 
the probability current density J in conjunction with the 
Schrodinger time evolution equation 



^ [^(x,y,t)ip(x,y,t)] +V- J = 0, 
d 



(6) 



where ^(x,y,t) is a two component wavefunction of the 
form tp(x,y,t) = (<f>i(y, t), <fci(y, t))e %k " x ] we obtain 



J x (y,t) = 2ARe[<t> 1 (y,t)(t>* 2 (y,t)]^ 

2Bk x [\Mv,t)\ 2 -\Mv,t)\ 2 ] 



(7) 



For k x — > 0, the first Dirac-like term dominates the sec- 
ond Schrodinger-likc term, implying that all low-energy 
edge states carry virtually the same current. The pro- 
file of J x is shown at the instants of time when the edge 
state exists on opposite edges in Fig. |i|b) . The sign of 
the current reverses when the state moves to the opposite 
edge. 

The evolution of the SHC carried by edge states of 
a given energy near a given edge (say, —L/2) can be 
obtained by integrating J x from —L/2 to —L/2 + 1/A. 
Since there are two oppositely propogating edge states 
at a given energy which exist on opposite edges of the 
system, we must add the currents carried by both of 
them. The time evolution of such a current following 
a sudden quench to the QCP is shown in Fig. [5] 
This current also displays a pronounced collapse and 
revival for sevaral cycles when the LE does so. The 
time-averaged value of the current is non-zero and is a 
significant fraction of the original value of the current. 
Thus, there is a persistence of the SHC carried by low- 
energy edge states following a sudden quench to the QCP. 

To summarize, the quantum fidelity shows a markedly 
different behavior for the 2-D DH in comparison to the 
1-D case. We observe interesting decoherence dynamics 
of the low-energy edge states of a TI Hamiltonian when 
quenched to the QCP displayed in the temporal evolution 
of their LE. Given the recent prospects of tuning of con- 
trol parameters of the related Hamiltonians, we believe 
that this study is of experimental relevance. 
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